# Figure S3: Alternative income threshold (correlation)

fig4_dat2 = bind_rows(import(here("data","fig4_inc2_agree_dat.rds")) |> 
                        mutate(agree = "Agree"),
                      import(here("data","fig4_inc_disagree2_dat.rds")) |> 
                        mutate(agree = "Disagree")) |> 
  invert_group()

fig_inc2_agree_scat = amce_scatterplot(fig4_dat2, "Income", "blue", -0.05, .45, "Agree")
fig_inc2_disagree_scat = amce_scatterplot(fig4_dat2, "Income", "blue", -0.05, .45, "Disagree")

fig_inc2_scat = fig_inc2_agree_scat + fig_inc2_disagree_scat

print(fig_inc2_scat)

ggsave(here("figures","appendix","figure_s3.png"), fig_inc2_scat,
       dpi = 600, units = "in", width = 10, height = 5)

ccat("See 'figure_s3.png' for Figure S3 output")